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We derive analytical expressions for the single mode quantum field state at the individual output 
ports of a beam splitter when a single-photon Fock state and a coherent state are incident on 
the input ports. The output states turn out to be a statistical mixture between a displaced Fock 
state and a coherent state. Consequently we are able to find an analytical expression for the 
corresponding Wigner function. Because of the generality of our calculations the obtained results 
are valid for all passive and lossless optical four port devices. We show further how the results can 
be adapted to the case of the Mach-Zehnder interferometer. In addition we consider the case for 
which the single-photon Fock state is replaced with a general input state: a coherent input state 
displaces each general quantum state at the output port of a beam splitter with the displacement 
parameter being the amplitude of the coherent state. 
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I. INTRODUCTION 

Quantum optics is an exciting field, in which many fundamental experiments, revealing the peculiarities of quan- 
tum mechanics, have been conducted. The advantage of optical experiments, compared with other fundamental 
experiments, is their simplicity. Quantum optics has been one of the main vehicles in the development of quantum 
information technologies and in particular of quantum cryptography [T]-[4] and optical quantum computing [3 [6] . 
The most widely used quantum states in this respect are coherent states and Fock states [7] . Important basic building 
blocks in quantum optics are beam splitters and Mach-Zehnder interferometers, which are passive and lossless four 
port devices [8-10 . Many cases of interference between vacuum, Fock states and coherent states have already been 
theoretically studied [10]. A famous example is the so called Hong-Ou-Mandel effect [11] where two single-photon 
Fock states arrive simultaneously at each input port of a balanced beam splitter. In this case there are no coincident 
photons at the output ports. On the other hand interference between a fluorescent photon and a classical field has 
been investigated [12 where the photon is created from a coherently excited atom going through its Rabi cycle of 
oscillation. This results in a time dependence of the interferometric fringe visibility as a function of the atomic Rabi 
frequency [T3] , 

Without going into details of the time-dependent photon generation a somewhat simpler but nevertheless important 
problem arises when we consider the interference between a true single-photon Fock |1) state and a coherent state \a). 
Based on an experiment in 2002 a displaced Fock state has been synthesized to a good approximation by overlapping 
a single-photon Fock state with a strong coherent pulse on a highly reflective beam splitter [14]. To the best of our 
knowledge an exact analytical expression for the evolving quantum state in this experiment has not been published 
so that the reconstructed state (Wigner function) of the beam splitter output could be compared with the theory. 

In this work we calculate analytically the output of a beam splitter (BS) in the case of a single-photon Fock state 
|1) and a coherent state \a) impinging respectively on its two input ports. In particular we derive an exact expression 
for the quantum state of the separate beam splitter output ports. Because of the generality of our calculations, the 
results we obtain are not restricted to the beam splitter, but are valid for all passive and lossless optical four port 
devices. In addition we consider the case for which the single-photon Fock state is replaced with a general input 
state. We show that a coherent input state displaces the quantum state at the output port of a beam splitter in phase 
space, whereby the displacement parameter is the amplitude of the coherent state. In the second part of this paper 
we show how the results can be adapted to the case of a Mach-Zehnder interferometer (MZI). Additionally, the mean 
photon numbers at the output ports of the MZI are determined. These results serve as a first step of investigating 
more complex systems that might find application in optical quantum computation. 

The intention of this article is to present the matter in a concise and didactic way. 
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FIG. 1. Beam Splitter (BS) with a single-photon Fock state |1) and a coherent state \a) 1 incident on its input ports. 



II. INTERFERENCE AT A BEAM SPLITTER 

We consider a general beam splitter with complex reflection coefficients r and r' and transmission coefficients t 
and t r . The phase relation between the coefficients depends on the construction of the beam splitter [15]. In the 
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Heisenberg picture the annihilation operators a of the incident fields transform as [10] 




where the indices indicate the corresponding ports (modes) . The unitary scattering matrix B must satisfy - for lossless 
devices - the so called reciprocity relations due to Stokes [16] 

\r'\ = \r\, \t\ = \t'\, \r\ 2 + \t\ 2 = 1 , r*t' + r't* =0. (2) 



A. Input: Fock state and Coherent state 

In our setup (see Fig.fTl) the incident field states are a single-photon Fock state |1) and a coherent state \a) t which 
can be written as the following product state: 

|1) la), = D^al |0> |0) x = e^""*^ |0) 10)! , (3) 

where is the creation operator and D\(a) is the unitary displacement operator. 
From Eq. ([I]) together with Eq. ([2| we easily obtain the following relations: 

aj = t'a\ + r' a\ and a[ = ra\ + ta\. (4) 
Obviously two vacuum states at the input ports of the beam splitter transform into vacuum states at the output 

B S 

ports: |0)q|0) 1 > |0) 2 |0) 3 . We use the Baker-Campbell-Hausdorff formula (TO] and Eq. Ml) to see how the input 

state Eq. ([3]) transforms into the corresponding output state under the action of the beam splitter: 

|!) I">1 e "(r4+t4)-^(r'&2+f& 3 )^^ 2 +r / a t) 

= D 2 (ra)D 3 (ta)(t'ai + r'a\) |0) 2 |0) 3 . (5) 
The density operator for the output state, Eq. ([5]), reads 

P23 = D2(ra)D 3 (ta)^ 23 Dl(ta)Dl(ra), (6) 

where 

p 2Z = (t'a\ + r'a\) |0) 2 |0) 3 (0| 3 (0| 2 (t'*a 2 + r"o 3 ). (7) 

In fact P23 would be the density operator of the beam splitter output with vacuum instead of the coherent input state. 
The two output modes (2 and 3) are in an entangled state. If we consider only output mode 3 we have to find the 
reduced density matrix by taking the partial trace over output 2, i.e. 



P3 



Tr 2 (p 23 ) = D 3 (ta) Tr 2 (^M^^M) I>t(ta) 
= D 3 (ta) Tr 2 (^(raJ^H^) r>t(ta) 

CO 

= D 3 (ta) ( n l 2 ^23 |n) 2 Dl(ta). (8) 

The operators D 3 (ta) and D\(ta) are not effected by the trace over mode 2, therefore we can put them outside the 
trace. Further, in the second and third line of the last equation we made use of the rule that the trace is invariant 
under cyclic permutations, Tr (ABC) = Tr(CAB), and the unity relation & (a)D(a) = 1. Inserting p 2 3 into the last 
equation we get with some basic boson algebra 

p 3 = \tf \ta) 3 (ta\ 3 + \r'\ 2 D 3 (ta) |1) 3 (1| 3 D 3 (to). (9) 

The result is a mixed state which is a convex combination of two pure identically displaced states: a coherent state 
(displaced vacuum) and a displaced Fock state [17] , where the displacement and the share of each pure state depends 
on the reflectivity and the transmittivity of the beam splitter. For the other output port we get the analogous result 

h = \r\ 2 \ra} 2 (ra\ 2 + \t\ 2 D 2 (ra) |1} 2 (1| 2 t)\(ra). (10) 
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B. Wigner function 



With the simple form of the density operator for the output mode 3, Eq. ([9|, we can calculate its Wigner function. 
The Wigner function is defined as [T8H2T] 



(ii) 



Wp(q,p) = J (q - y/2\ p\q + y/2) exp(iyp/h)dy, 



where p and q are the field quadratures. By inserting Eq. (|9| into Eq. (11) we see immediately that we get a sum of 
two Wigner functions at the output port, 



W h = \t'\ 2 W Kb{tam) + \r'\ 2 W mtamy 



(12) 



where W 



p(D(ta)\0)) 



is the Wigner function of the coherent state \ta) and W^^^^ is the Wigner function of the 



displaced Fock state D(ta) |1). These two individual functions are well known [22] [17]. For the coherent state it is 
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W: 
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p(D(ta)\o))(^P) = W pm) (q,p) - ^ 



exp 



PQo 

h 



(13) 



and for the displaced Fock state it is 



WD(ta)\l))(^P)= W K\l))(<liP) 



exp 



PQo 

h 



1 + 2 




where qo = y/h/uj. The displacement quadratures 



q' = q — V2qo(Re ta) and 



p — p 



a/2 — (Im ta) 
Qo 



(14) 



(15) 



shift the minimum of the Fock state W, 



-1 (and the maximum of the coherent state W/s(|o)) = 1) to g 



>/2go(Re ta) and p = \/2^(Im ta). 

The Wigner function Eq. ( [12] ) for the output field at port 3 is depicted in Fig. [2] for two different cases. In the first 
case we consider a 50:50 BS, where the Wigner function is an equal mix between a coherent and a displaced Fock 
state. The absolute minimum of this function is zero, situated at the center of the displaced Fock state, as can easily 
be seen from Eqs. (12)-(14), and the total Wigner function is non-negative. In the second case we consider a highly 
reflective (99:1) beam splitter comparable to the experiment mentioned in the introduction [H] , in which a displaced 
Fock state has been synthesized. However the high reflectivity would lead to a smaller displacement at port 3, unless 
the intensity of the incoming coherent state is increased to get an identical factor ta for Figj2^i and Figj2}}. Therefore 
the shift is identical, but due to the different reflection and transmission coefficients of the BS the relative weights of 
both input states are different. Indeed we see from the Figj2}} that the displaced Fock state is very dominant with a 
pronounced minimum nearly reaching the original value of -1. 

We have thus established an exact theory for the experiment of creating a displaced Fock state [14] instead of the 
theory which the experiment was originally based on [23-25 or the theory in [26] . 



C. Input: Coherent state and a general state 



For future theoretical considerations and experimental demonstrations it might be useful to replace the ideal Fock 
state by a general state 



The input state then reads 
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FIG. 2. Wigner function (X7r, h = uj — 1) of the beam splitter output 3, Eq. (12). It is a statistical mixture between a coherent 
state and a displaced Fock state, (a) 50:50 beam splitter (t = l/y/2) with the coherent input state \a = v^e 27 ^ 4 ). (b) Highly 
reflective beam splitter (t = 1/10) with the coherent input state \a = 10e 27r/ ^ 4 ). Note, that ta is identical for both (a) and (b). 



It is straightforward to see that the output mode 3 of the BS is given by an equation analogous to Eq. ([8|, whereby 
the density matrix p23 becomes 

oo oo 1 

= E E -r^c m ct(t'al + r'atr |0} 2 |0) 3 (0| 3 (0| 2 (t'*a 2 + r'*a 3 ) 1 . (18) 

Similar to above P23 is the non displaced density operator which would emerge at the output if the coherent state at 
the input would be replaced by vacuum. Thus the coherent input state displaces any state at the individual output 
port of a beam splitter compared to a vacuum input (see Eq. In the limit of a highly reflective beam splitter 



(t — » 0) Eq. (18) becomes /?23 — >• |0) 2 1^)3 {^3 (0 1 2 . If we furthermore consider a strong coherent state (ta finite), we 
see from Eq. that the effect on an arbitrary input state is approximately a mere displacement by ta of this state. 
A similar result has already been shown using a different approach in [26] confirming our calculations. 



III. INTERFERENCE AT A MACH-ZEHNDER INTERFEROMETER 



Now we consider a single-photon Fock state and a coherent state at the input of a Mach-Zehnder interferometer 
(MZI) as shown in Fig. [31 Again we want to calculate the quantum state at the output ports. 



A. MZI scattering matrix 

As suggested in [9], a MZI can be considered as a lossless and passive four port device. The annihilation operators 
of the field states transform similar to Eq. ([T]) and yield 

ft)= B -'(":)- < i9 > 

The necessary and sufficient condition for the scattering matrix Bmzi is that it has to be unitary. In fact, since 
we have already considered a general scattering matrix in the case of the beam splitter, Eq. ([T]), which essentially 
accounts for all passive and lossless four port devices, the calculations from the last section are generally valid. The 
scattering matrix of an arbitrary four port MZI, in particular, may be represented as a composition of several unitary 
scattering matrices. 
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FIG. 3. Mach-Zehnder Interferometer (MZI) with a phase shift in path 3. The input states at BSi are a single photon Fock 
state |1) and a coherent state \a) 1 . 



The scattering matrix for the MZI, depicted on Fig (3] is composed of 3 unitary matrices, two matrices corresponding 
to the two beam splitters (see Eq. 0) and one matrix associated with the phase shift, which in our case reads 



U e 



1 ( 

e 



Accordingly the annihilation operators of the field states transform as 
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where in the last line we defined the following variables: 



tm = nt 2 + e l tir' 2 , 
tM = rir 2 + e l6 tit' 2 , 



' M 



1'2- 



(20) 



(21) 



(22) 



The variables defined in Eq. (22) obey the reciprocity relations Eq. ([2]), in particular |t^| = \t' M \ and \tm\ 



'Ml 



Here we add a short comment with respect to Figj3j The output ports 4 and 5 can serve as two paths of a second 
MZI including a second phase shift, subsequently followed by a third MZI of a similar type. Such a cascade of 3 
interferometers can operate as preparation, distribution and measuring module of a simple linear optical gate helpful 
in quantum computing systems. The resulting scattering matrix of such a device would be a matrix product of 3 
matrices similar to Bmzi used in Eq. (21). 



B. MZI Output states 

The states at the output ports 4 and 5 of the MZI can now be calculated analogously to the output states of 
the single beam splitter in the previous section. Therefore according to Eq. ^ and by replacing the beam splitter 
scattering matrix B with the MZI scattering matrix Bmzi we get for the MZI output state 



!)o D A (r M ^)D b (t M (y)it' M a\ + r' M a\) |0) 4 |0) 



(23) 
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Further we get for the reduced density matrix at output port 5 (compare with Eq. (I9|), 



and output port 4 



h I^m| 2 \tMOL) h (t M a\ 5 + Vm\ 2 b$(t M ot) |1> 5 (1| 5 bl(t M a) 



Pa = \vm\ \rMOt) 4 {tmol\ a + \t M \ D A {r M ot) |1) 4 (1| 4 £>4(r M a)- 



(24) 



(25) 



The result is again, as in the case of the beam splitter, a mixed state between a coherent state (displaced vacuum) 
and a displaced Fock state, both displacements being identical. The Wigner function of this state can be calculated 
using the analogue of Eq. (12). The calculations where a coherent state and a general state are incident on a MZI are 



equivalent to the beam splitter case, so the output state is again merely displaced compared to a MZI with vacuum 
input instead of the coherent input state. 



C. Photon numbers 

Finally we calculate the average photon number at the output ports of the MZI. For port 5 it is defined as 



t M OL) {t M ® 



n=0 



n ) = {tMOL 



(n) 5 = Tr 5 (n5/95) = ^ (n\ a\a*>fe 
For the coherent state we obtain, of course 

oo 
n=0 

For the displaced Fock state we obtain 

CO 

^(n\a\a b D b {t M a) I) (l b\{t M a) 



n 



tMOL) = \t M ot\ . 
5 



(26) 



(27) 



= (l\Dl(tMa)ala 5 D 5 (t M a) l) =l + \t M a 



(28) 

where in the last step the commutation relation [7 aD(tM&) — D(tMOL){a + tM<x) nas been used. With the last two 

1 2 . i 1 2 



equations and using that |£m| + \ t m\ = 1 we can now calculate the average photon number, Eq. (26), and get 

(n} 5 \t M \ 2 \tMOL\ 2 + \r M \ 2 (l + \tMOi\ 2 ^j = \r M \ 2 + |^m«| 2 • 
For port 4 we get analogously 



(29) 



(n) 4 = \t M \ + \rM<x\ ■ (30) 
Using this formalism the mean square deviation ((An) 2 ) = (n 2 ) — (n) 2 of photons can be calculated leading to 

((An) 2 ) 5 = \t M a\ 2 (l + 2|r M | 2 ) + \r M \ 2 \t M \ 2 , 

((An) 2 ) 4 = |r M a| 2 (l + 2\t M \ 2 ) + \r M \ 2 \t M \ 2 (31) 

for the two output ports. 

D. Discussion of a MZI with balanced BSs 



As a particular example, that can be realized easily in an experiment, we calculate the average photon numbers 
for the case of two 50:50 dielectric layer beam splitters with a phase factor of e t7V / 2 = i for the reflected beams 



(h = t[ =t 2 = t' 2 = l/y/2, 



7*1 



r 2 



i/V2) [E]. In this case the transmittivity and the reflectivity are 



|^m| 2 = sin 2 (f9/2) and \tm\ 2 = cos 2 (f9/2). Both are completely determined by 0. The average photon numbers are 
then 



(n) 4 = sin 2 (<9/2) + \a\ 2 cos 2 (<9/2), 



(32) 
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(n> 5 = cos 2 (0/2) + \a\ 2 sin 2 (0/2) (33) 
and the mean square deviations are given by 

<(An) 2 > 4 = \ sin 2 (61) + M 2 [cos 2 (^) + \ sin 2 (6>)] , ((An) 2 ) 4 - ((An) 2 ) 5 = |«| 2 cos(0) . (34) 

Figure [4] shows the average photon numbers as a function of the phase shift for different coherent input states \a). 
The maxima and minima in the graph indicate the cases where the quantum states at the separate output ports are 
either a pure coherent state or a pure Fock state. 

<n> 4 <n> 5 
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FIG. 4. Average photon number at output 4 (a) and output 5 (b) of the Mach Zehnder Interferometer with the coherent input 
states 1 1 a | = 0), ||a| = 1) and ||a| = 2). Note that for \a\ 2 — 1 the mean photon numbers are independent of and equal to 1. 

The technical advantage of the MZI (with two 50:50 BSs) over a single beam splitter is, that tM and vm can simply 
be adjusted by changing the phase in path 3. All desired states could be generated in experiments without changing 
the setup by selecting another beam splitter with the appropriate transmittivity |t| 2 and reflectivity |r| 2 . 

Let us compare the result for \a\ 2 = 1 with the case when two separate single-photon Fock states are impinging on 
the input ports of an MZI. The first balanced BS± of the latter is reminiscent to the Hong-Ou-Mandel (HOM) effect 
[TT] . After the BS\ the state is z ( 1 2) 2 1 0)3 + |0)2 |2)s) / a/2 which means that no coincidences appear behind BS\ [TO] , 
Introducing a phase shift 6 in path 3 and a second BS2 (see fig. [3]) the wave function behind the MZI reads 

{(1 - e™)[ |0> 4 |2>5 - |2> 4 |0> 5 ] + V2(l + e 2 ^)|l) 4 |l) 5 } • (35) 



2y/2 



The probability for coincidences is W45 = cos 2 (6) while the probabilities of measuring 2 photons at the output ports 
4 and 5 are W4 = W5 = \ sin 2 (0). However, the mean photon numbers (n)^ and (n)^ are equal to 1 and therefore 



independent of 0. This can be compared with the results of the present paper (see Eqs. (32) and (33)) in which the 
mean photon numbers in output 4 and 5 also do not depend on 0. Anyway, one has to keep in mind that these two 
examples have completely different initial conditions although the input mean photon numbers are the same in case 



of \a\ 2 = 1. In fact, one has to compare the entangled quantum states of the HOM-effect Eq. (35) on the one hand 



with the MZI-state Eq. (23) on the other hand and realize that these are completely different. By executing the 
trace operation for the separate output ports, however, entanglement is eliminated in both cases yielding the identical 
outcomes for the mean photon numbers. 



IV. SUMMARY AND OUTLOOK 



To sum up, we have derived simple analytical solutions for a quantum state and its Wigner function at the output 
ports of a general passive and lossless optical four port like a beam splitter, when a coherent state and a single-photon 
Fock state are incident on the input ports. These calculations could be of interest in the fields of quantum cryptography 
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based on continuous variables and optical quantum computing. We have obtained a statistical mixture between of a 
coherent state and a displaced Fock state and have derived the corresponding Wigner functions. Furthermore, we have 
shown that a coherent input state displaces the quantum state at the output port of a beam splitter in phase space 
as compared with vacuum at the input port, when in both cases a general state is incident on the second input port. 
Additionally, we have analyzed the quantum states behind a Mach-Zehnder interferometer and have evaluated mean 
photon numbers at the output ports. It turns out that for an input state |l)o|M = l)i the mean photon numbers 
behind the MZI do not depend on the phase shift 9 inserted in the interferometer. 

In addition to the 3-interferometer-cascade mentioned in section III/ A a further possible application of the presented 
formalism can be proposed. Instead of the reflector R in path 2 of Fig. [3] a beam splitter can be inserted in such a 
way that beam 2 is both reflected and transmitted. The beam splitter provides also an additional input. A second 
MZI then can be added below the first one. In particular these two interferometers have one beam path in common 
generating a so-called two-loop interferometer which consists of 4 beam splitters and 2 reflectors providing 3 input and 
3 output ports. Such a device enables the superposition of three wave functions and will be theoretically investigated 
in a next step. 

The results of this paper can serve as a basis for further phase-space investigations of higher Fock states in com- 
bination with coherent, squeezed or thermal states as input states used in current optical setups. Those setups (BS, 
MZI, Phase gate, CNOT gate) are building blocks for linear optical quantum computers. 
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